AD-A102  342 


UNCLASSIFIED 


AERONAUTICAL  RESEARCH  ASSOCIATES  OF  PRINCETON  INC  NJ  F/G  20/4 

STATISTICAL  CONSTRAINTS  ON  SCALAR  VARIABLES  IN  TURBULENT  FLOWS* (U) 
FEB  81  G  SANDRI#  P  J  MANSFIELD »  A  K  VARMA  F44620-76-C-0048 

ARAP-346  AF0SR-TR-81-0593  NL 


AFOSR-TR-  81  -0598 


V8  j 

fW 


A.R.A.P.  Report  No.  346 


STATISTICAL  CONSTRAINTS  ON  SCALAR  VARIABLES 
IN  TURBULENT  FLOWS 


Aeronautical  Research  Associates  of  Princeton,  Inc. 
50  Washington  Road 
P.0.  Box  2229 

Princeton,  New  Jersey  08540 


February  1981 


Scientific  Report 


The  views  and  conclusions  contained  in  this  document  are 
those  of  the  authors  and  should  not  be  Interpreted  as 
necessarily  representing  the  official  policies  or  endorse¬ 
ments,  either  expressed  or  Implied,  of  the  Air  Force  Office 
of  Scientific  Research  or  the  U.S.  Government. 

JLpprreed  tv*  pMi*  vela 

atntMbutfon 

Prepared  for 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
Building  410 

Bolling  Air  Force  Base,  DC  20332 


Si  b 


Oi> 


*0 


Unclassified _ 

security  CLASSIFICATION  or  This  PACE  fWhmn  Da  10  Enfitd) 

f ’ " .  REPORT  DOCUMENTATION  PAGE 


2.  GOVT  ACCESSION  NO. 


It l  AFCfSR.TR^8  1  -0593  4x/li 

«.  Title  cmd  Submit) 

j  Statistical  Constraints  on  Scalar  Variables  4 
fa  in  Turbulent  Flows,  [./ 


7  AUTHOR/*)  ^ 

Gui do /Sandri/  Peter  J.^Mansfield^"  V 

10,  Ashok'  K . yvarma ~gprir  Co  reman  duP.y  Donaldson  ' ' 

9.  performing  organization  name  and  address 

Aeronautical  Research  Associates  of  Princeton, Inc. 
50  Washington  Road,  P.0.  Box  2229  // 

Princeton,  New  Jersey  08540 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Research /NA  f //  . 
Bolling  Air  Force  Base,  DC  20332 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

J  RECIPIENT'S  C  ATALCkTnuMSER  ~ 


5.  type  or  report  a  period  covered 

Scientific  ^  ^ 

1  Jan  7R.  — •  31  Jul  78  v  _  _ 

I  FeftPORMNO  ORG.  REPORT  WOm»CR 

A.R.A.P.  Report  No.  346 

i  CONTRACT  OR  GRANT  NUMBER!.; 


>^446^0-76- C-, 


to.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  *  WORK  UNIT  NUMBERS 


2307A2  / 
6T102F 


/l) 


February  1981 


14.  MONITORING  AGENCY  NAME  A  ADDRESS (U  dUtmtmxt  from  Controlling  Oiiico)  |  15.  SECURITY  CLASS,  (of  */»/•  "port) 


)  )  r  P- 


[  16.  DISTRIBUTION  STATEMENT  (of  thio  Report) 


Unclassified 


IS*.  DECLASSIFICATION 'DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  tho  abstract  entered  In  Block  20,  If  <Wf/«renl  from  Roport)  £  *T  | 

V' 


x- 


[  10.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  ravarta  aide  If  neceaemry  and  Identify  by  block  number) 

Turbulence 

Statistical  fluctuations 
Mixing  layer 

Turbulent  chemical  reactions 


20.  ABSTRACT  (Continue  on  rovoroo  oldo  If  neceaemry  and  Idontlfy  by  block  number) 

*  We  consider  the  statistical  behavior  of  scalar  variables  in  turbulent  flows. 
The  presence  of  chemical  reactions  requires  statistical  models  for  third  and 
higher  moments  in  order  to  close  the  rate  equations.  A  realistic  analysis 
should  not  be  restricted  to  small  fluctuations.  The  results  presented  here 
are  completely  free  from  such  restrictions.  The  essential  points  to  be 
demonstrated  are: 

1)  Given  'A)  and  *  ,  we  can  obtain  stringent  statistical  bounds  on  A1 zBr 
and  related  moments.  These  bounds  are  found  to  be  of  interest  in  discussing 


DD  ,55 


_ Unclassified _ O l?  £_j/  O  (0* 

SECURITY  CLASSIFICATION  OF  This  PAGE  om  Cm* ttH)  .  / 


.  ■  -  fc/j 


security  classification  or  this  pagc(T»Ti«o  d«<«  Swiwq 

recent  experiments. 

2)  Maximum  and  minimum  third  (and  higher)  moments  can  be  reached  only  with 

Dirac  functions  (i.e.,  discrete  distributions).  /•  •  *  ' 

3)  We  can  always  realize  a  statistically  acceptable  choice  of  “K  and  A'2 
with  a  few  Dirac  functions  (the  minimum  number  is  two  distinct  ones). 

The  statements  given  above  hold  true  when  several  moments  (rather  than  only 
and  Ari  are  given  as  well  as  when  the  means  of  several  variables  are 
given.  As  one  consequence,  we  shall  see  that  we  can  find,  for  the  purposes 
of  modeling,  a  discrete  distribution  that  represents  the  desired  set  of 
moments  and  that  is  statistically  legitimate  for  all  allowed  values  of  the 
derived  moments  (box  model). 

s 


\  CIIC  lies  ° _ 

\  V*'"**'^  Co4'3. 

V  d/or 


1 


Unclassified 


SECURITY  CLASSIFICATION  OF  This  PAOC<m»n  Dm*  Cnl*r*tf; 


TABLE  OF  CONTENTS 


I.  Introduction  1 

II.  Statistical  Analysis  5 

III.  Incompressible  Flow  with  Variable  Density  17 

IV.  General  Method  for  the  Determination  of 

Statistical  Bounds  23 

V.  Model  Selection  37 

Appendices 

1.  Convexity  of  Moment  Domains  45 

2.  Standard  PDF's  47 

3.  Relations  Imp! ied  by  the  Equation  of  State  51 

4.  Application  of  Instantaneous  Bounds  to 

Ensure  Correct  Statistics  for  Solutions 

of  the  Rate  Equations  53 

Figures  1  through  50  61  ff. 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  fAFSCl 

NOTICE  OF  TRANSMITTAL  TO  DDC 

This  technical  report  has  boon  roviewed  and  is 

PUtlic  rk;lea5e  I  AW  AFR  190-12  (7b). 
Distribution  is  unlimited, 

A.  D.  BLCSE 

Te jhnl Ofi  1  In'ormation  OiTicer 


I.  INTRODUCTION 


We  consider  below  the  statistical  behavior  of  scalar  variables  in 
turbulent  flows.  The  presence  of  chemical  reactions  requires  statistical 
models  for  third  and  higher  moments  (e.g.,  A'2B'  and  K ' A '  B‘)  in  order  to 
close  the  rate  equations.  A  realistic  analysis  should  not  be  restricted  to 
small  fluctuations.  The  results  presented  here  are  completely  free  from 
such  restrictions.  The  essential  points  to  be  demonstrated  are 

— >r 

_ 1)  Given  A  and  A'  ,  we  can  obtain  stringent  statistical  bounds  on 

2 

A'  B'  and  related  moments.  These  bounds  are  found  to  be  of  interest  in 
discussing  the  CalTech  experiments  (e.g.,  as  reduced  by  Konrad*). 

2)  Maximum  and  minimum  third  (and  higher)  moments  can  be  reached  only 
with  Dirac  functions  (i.e.,  discrete  distributions). 

3)  We  can  always  realize  a  statistically  acceptable  choice  of  A  and 

2" 

A'  with  a  few  Oirac  functions  (the  minimum  number  is  two  distinct  ones). 

The  statements  given  above  hold  true  when  several  moments  (rather  than 
-  1 

only  A  and  A1  )  are  qiven  as  well  as  when  the  means  of  several  variables 
are  given  (e.g.,  A,  A  ,  a  ,  o'^  ).  As  one  consequence,  we  shall  see  that 
we  can  find,  for  the  purposes  of  modeling,  a  discrete  distribution  that 
represents  the  desired  set  of  moments  and  that  is  statistically  legitimate 
for  all  allowed  values  of  the  derived  moments  (box  model). 

The  role  of  the  higher  moments  can  be  illustrated  with  a  second-order 
chemical  reaction 

A  +  B  C 

We  have  the  familiar  rate  equation 

ar  *  ar  -  -k<t)ab 

It  follows  that  the  product,  AB  ,  satisfies 

(AB)  *  -K(AB2  +  A2B) 

In  terms  of  the  Reynolds  decomposition 

A  -  A  +  A' 

we  find,  in  the  absence  of  temperature  fluctuations, 

♦Konrad,  J.H.,  "An  Experimental  Investigation  of  Mixing  in  Two-Dimensional 
Turbulent  Shear  Flows  with  Applications  to  Diffusion-Limited  Chemical 
Reactions."  Graduate  Aero. Labs.,  California  Institute  of  Technology,  Pro¬ 
ject  SqUID  Tech.  Report  CIT-8-PU  (1976) 


^  (AB  +  A ' B ' )  =  -K(A?  +  S) 

where 

A?  =  <^(A  +  A')(B  +  B')2^> 

=  A  B2  +  A  ?  +  2B  Fb  +  A'B'2 

— 7 

We  shall  use  below  this  third  moment,  AB  ,  to  test  the  stringency  of  the 
statistical  bounds  derived. 

When  temperature  fluctuations  are  included,  a  joint  distribution  in  T, 

A  and  B  need  be  assumed  to  close  the  rate  equations.  Although  we  do  not 
develop  temperature  models  in  detail  in  this  report,  we  note  that  many  of  the 
statistical  (lower)  bounds  derived  for  A  do  apply  to  the  temperature 
fluctuations;  for  example,  the  third  moment  of  the  concentration  satisfies 

A^  >  ^  (A^  -  A2) 


for  any  0  <  A  £  1  and  0  <  A’^  <  A(1  -  A);  similarly,  the  temperature 
third  moment  satisfies 

— ^  2 


>  L_  (r7  -  t2) 

T 

— 7  — 3 

for  any  T  and  any  V  positive.  The  usefulness  of  the  A'  bound  is 

shown  below  by  ruling  out  the  model  A'  =0  over  the  edges  of  a  mixing 

layer. 

The  statistically  allowed  values  of  anj/  set  of  moments  can  be  shown 
(Appendix  1)  to  form  a  domain  that  is  necessarily  convex;  therefore,  the 
domain  has  no  "holes."  This  property  is  very  useful  when  a  numerical  scan 
is  performed. 

The  statistical  bounds  that  we  obtain  are  "sharp”;  that  is,  we  can 
prove  in  each  case  that  the  equal  sign  is  realized  by  a  (unique)  probability 
distribution  that  is  discrete.  It  follows  that  the  bounds  obtained  are 
"best  bounds"  and  cannot  be  improved. 

Although  we  apply  here  our  analysis  to  mixtures  with  two  species  (He 
and  Ng).  the  results  derived  for  the  (molar)  concentration  A  (and  its 
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conjugate  B  =  1  -  A)  apply  to  each  concentration  of  an  n-component 
mixture. 

We  develop  a  model  for  the  probability  of  the  molar  fraction  by 
choosing  a  discrete  distribution  that  realizes  the  mean  and  second  moment  of 
the  molar  fraction  as  well  as  the  mean  and  second  moments  of  the  mass 
fraction.  The  model  selected  contains  a  free  parameter.  Since  the  empir¬ 
ical  evidence  is  that  the  mean  entropy  of  mixing  is  below  the  maximum 
allowed  statistically  (in  qualitative  agreement  with  the  dominance  of  large 
scale  structures  in  the  turbulent  mixing),  we  fix  our  free  parameter  so  as 
to  minimize  the  mean  entropy  of  mixing. 


3 


II.  STATISTICAL  ANALYSIS 


To  illustrate  the  main  features  of  our  statistical  analysis,  we  discuss 
below  in  detail  the  behavior  of  the  second  moments  and  we  give  results  for 
the  third  and  fourth  moments. 

We  introduce  the  analysis  with  a  useful  diagram,  the  "box"  diagram  which 
gives  a  concise  representation  of  a  (one-dimensional)  probability 
distribution. 

Given  a  P(A)  >  0  ,  the  quantity 


t(Ao) 


P(A)  dA 


is  the  fraction  of  time  during  which  A  <  Aq  .  The  graph  of  the  inverse  of 
the  map  x(A)  is  the  box  diagram.  For  example, 

P(A)  =  A6(A  -  1)  +  (1  -  A)6 (A)  (2 

gives  x(A)  =  (1  -A)H(A)  +  AH(1  -  A)  where  the  heaviside  function  is 


H(A)  * 


1  A  >  0 
0  A  <  0 


and  then. 


0  0  <  t  <  1  -  A 
1  1  -  A  <  x  <  1 


(see  Fig.  1). 

We  will  now  consider  the  bounds  on  the  second  moment  for  a  given  value 
of  A  .  We  first  observe  that 


0  <  A  <  1 


implies  the  two  inequalities 


0  <  AP(A)  ;  AP(A)  <  P(A) 

because  P(A)  >  0  .  Integrating  both  between  0  and  1, 


0  <  A  <  1 


HiECBttLNQ  na  KJUK-NOT  FIJuMH) 


By  a  similar  reasoning,  we  can  prove 

A2  <  /?  <  A 

To  prove  the  upper  bound,  we  have  successively 

by  (4U),  i.e.,  the  upper  bound 
by  (4L),  i.e.,  the  lower  bound 
because  P  >  0 ,  q.e.d. 

To  prove  the  lower  bounds,  we  use  Reynolds'  decomposition 

A  =  A  +  A' 


A  <  1 
A2  <  A 
A2  <  A 


(7) 


(8) 


(9) 


Thus , 


?  ■  <^( A  +  A ‘  )2y  =  A2  +  A*7  >  A2  because  A'2  >  0,  q.e.d. 


The  bounds  on  the  second  moment  have  been  obtained  by  elementary  methods. 
We  can  similarly  show  that  the  bounds  (A2"  =  A  and  A2  =  A2)  are  realized 
by  unique  discrete  distribution. 

For  the  upper  bound,  we  have  successively 

~7 

A  =  A 


and 


/. 


A  -  A  =  0 


1 


(A  -  A  jP(A)  dA  =  0 


However,  the  integrand  is  positive  (or  zero)  by  Eq.  (8)  and  hence  the  inte¬ 
grand  must  vanish: 


(A  -  AC)P( A)  =  0 


(10) 


If  P  is  to  be  nonzero,  we  must  have  A  =  0  or  A  =  1  and,  therefore. 


P(A)  =  (1  -  e)fi(A)  +  e5(A-  1) 


(ID 


*> 
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Since  A  is  given,  e  is  determined.  In  fact,  by  direct  computation, 
A  =  e  and 


P( A)  =  (1  -  A)6(A)  +  A6(A  -  1 ) 


(12) 


which  is  the  distribution  given  in  Eq.  (2)  above.  We  note  that  the  P(A) 
that  realizes  the  upper  bound  of  A^  ,  Eq.  (12), is  uniquely  determined  by  A  . 
We  can  view  relation  (12)  as  a  consequence  of  (10)  and  of  L.  Schwartz1 
theorem  x6(A)  =  0  . 

The  lower  bound  is  also  realized  by  a  unique  distribution.  We  have 
from  (7L) 


or 


-  A2)P(A)  dA  = 


(A  -  A)2P(A)  dA  = 


1 

A,2P(A)  dA 


As  for  the  upper  bound,  the  integrand  is  positive  or  zero  and,  therefore, 

A,2P(A)  =  0 
or 

(A  -  A)2P(A)  =  0 
It  follows  that  A  =  A  if  P  ^  0  ,  or 


P(A)  =  <5{A  -  A) 


(13) 


This  distribution  can  be  represented  by  a  box  diagram  as 


<5(A  -  A)  = 


~7 

This  result  completes  the  proof  that  both  extreme  values  of  A  are 
realized,  for  given  A  ,  by  unique,  discrete,  probability  densities. 

A  summary  of  the  relations  obtained  so  far  is  obtained  in  describing 
the  statically  allowed  values  of  A  and  7  as  a  domain  in  a  plane  with 


(14) 
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-  T  -  ~2 

coordinates  A,  A  .  The  domain  of  moments  M^(A,  A  )  is  shown  in  Figure 

2  that  represents 


.{o<  A  <  1  ;  A2  <_7<_  a] 


We  note  that  the  maximum  width  of  the  domain  occurs  at  A  =  1/2  and  is 
W12  =  ^  =  •  Clearly  the  maximum  width  of  the  statistically  allowed 

domain  should  be  large  compared  with  the  experimental  error  if  a  statisti¬ 
cally  significant  experiment  is  contemplated.  We  show  below  that  the  Cal 
Tech  experiments  amply  fulfill  this  requirement. 

7 

We  recast  our  result  in  terms  of  the  centered  moment  A'  by  means  of 
the  identity 

— 7  ~7  -? 

A'*  -  hc  -  (v  (16) 

-  — 7 

The  domain  of  statistically  allowed  (A,  A'  )  pairs  is  readily  shown  to  be, 
from  Eq.  (15) 

Mi2  =  {°  -  ~k  -  1  ;  0  -  M1  '  A)}  (17) 

In  Figure  3,  we  show  a  diagram  of  .  We  illustrate  our  theoretical 

bounds  with  the  Konrad  data  in  Figure  4.  In  this  and  similar  following 
figures  of  this  section 

B  =  concentration  (18) 

— ^ —  =  similarity  variable  across  the  flow  (19) 

x  "  o 

A  =  1  -  =  0  (The  He  is  premixed  with  A  to  (20) 

rc2  achieve  molecular  mass  equality) 

s  =  1  =  Konrad's  mass  ratio  parameter 
r  =  .38  =  velocity  ratio  (He  is  faster) 


In  view  of  the  fact  that  B  is  monotonic  in  the  similarity  variable 
y/(x  -  xQ)  ,  we  may  plot  the  experimental  data  on  the  diagram  of  Figure  2 
using  the  B  axis  as  a  "stretched"  similarity  variable.  The  resulting 
diagrams.  Figures  4  and  5,  show  two  alternative  (and  equivalent)  representa¬ 
tions  of  the  second  moment  of  the  N2  concentration  &  and  of  its 
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statistical  bounds  as  dictated  by  the  corresponding  values  of  B  .  Konrad 

gives  pluts  of  both  the  Gaussian  measure  for  the  second  moment  J B ' ^  and  of 
the  (second  order)  mixedness  parameter 

l/VB7  1  _  A*7  ,  7  -  A2  (21) 

A  B  A( 1  -  A)  A(1  -  A) 

Readings  of  the  reproducibil ity  of  the  similarity  profiles  and  comparison 

2 

of  values  obtained  for  the  moment  B  and  B  from  the  distribution 
functions  vs  those  obtained  from  the  Konrad  plots  indicate  an  "experi¬ 
mental  error  bar"  (or  uncertainty)  of  ±0.01  on  the  scales  indicated. 
Sensitivity  of  the  bounds  to  errors  in  B  is  estimated  by  setting  B  = 

B  ±  0.01  .  Then  at  y/(x  -  xQ)  =  0  (center  of  flow), 


6L2  =  ±0.011 

(22) 

6U2  =  ±0.010 

(23) 

6W2  =  +0.001 

(24) 

where 

W2  =  U2  -  L2  <  0.25 

(25) 

We  see  that  the  width  of  the  statistical  bounds  is  less  sensitive  to  error 
than  the  lower  and  upper  bounds  themselves.  The  equivalent  shift  in  the 
similarity  variable  near  y/(x  -  x  )  *  0  is  estimated  to  be 
<5(y/(x  -  xQ))  *  -0.004. 

In  conclusion,  statistical  analysis  of  the  second  moment  indicates  that 
the  experiment  is  clearly  statistically  significant  and  that  the  second- 
order  fluctuations  are  on  the  low  side  in  the  central  portions  of  the  flow. 

As  a  final  point  in  the  statistical  analysis  of  the  second  moment,  we 
give  a  constructive  procedure  whereby,  given  any  allowed  value  of  the  pair 
(A,  A7)  (or  alternatively  of  (A,  A'  )),  a  discrete  distribution  is  found 
that  realizes  the  given  values.  The  construction  is  based  on  the  results 
obtained  on  the  realizability  of  the  extreme  values  of  the  moments.  In 
particular,  we  have  demonstrated  that  maximal  fluctuations  are  achieved  by 
discrete  distributions.  Exploiting  the  convexity  of  the  domain  of  values 
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- 

of  (A,  A  )  proven  in  App.  1,  we  proceed  as  follows  (Fig.  6).  At  points 
a  and  b  we  must  have 


P,  =  £<S(A)  +  (1  -  e)6(A  -  1) 

a 


(26) 


Pb  =  6(A  -  Aq) 

We  now  form  the  (convex)  combination 

P  •  lPa  ♦  (1  -  *)Pb 


where 


0  <  X  <  1 


(27) 

(28) 
(29) 


The  resulting  distribution  is  given  by  the  box  diagram,  Figure  7.  Straight¬ 
forward  calculation  (best  carried  out  with  the  general  method  of  the  char¬ 
acteristic  equation  given  in  Section  IV)  yields 

A2  -  AA„ 


1  -  A, 


e2  =  1  -  A  - 


A  -  7 


(30) 


A  -  A^” 

e3  =  A3(l  -  A3) 


where  A3  is  bounded  above  and  below  from  the  requirements,  respecti vely, 
that  Ej  be  positive  and  e2  be  positive 


B  A 


(31) 


Three  special  cases  are  particularly  interesting  since  they  correspond 
respectively  to:  (1)  maximal  A  moments;  (2)  minimal  A  moments;  and 
(3)  symmetric  A  and  B  moments.  The  three  special  cases  are  given  in 
full  in  Figure  8.  The  corresponding  cells  are  obtained  from  Figure  6. 

To  discuss  moments  higher  than  the  second,  we  need  some  additional 
tools.  These  are  incorporated  in  the  following  inequalities: 

1.  Holder's  inequality,  f  and  g  are  arbitrary  and  P  is  positive. 


Then, 


T 


ll/P 


fP(a)P(a)  da 


r 


il/q 


gq(a)P(a)  da 


/' 


provided 


or 


1/P  +  1/q  =  1 


7p  '/P  -q  ‘/I  _ 

f  •  gH  >  Tg 


1  |  ^(oi)g(ot)P(ot)  da 
0 

(32) 


(33) 

(34) 


Equality  holds  if,  and  only  if. 


fp  =  X  gq 


(35) 


Furthermore,  the  Cauchy-Schwa rz  inequality  corresponds  to  p  =  q  =  2;  that 
is 


with  equality  at 
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f  =  Xg 


(36) 

(37) 


The  conditions  for  equality  are  essential  in  our  analysis  since  they  allow 
us  to  determine  the  realizations  of  the  extreme  values  of  the  moments. 

2.  Jensen's  inequality.  If  ^(a)  is  convex  ( i . e . ,  <p"  >  0  throughout 
the  domain  of  interest  (0  <a  <_  1)),  then 


4>(a)  >_  4>(a) 


(38) 


3.  Tchebytcheff 's  inequality.  If  f  and  g  are  similar  (both  increase 
or  both  decrease  in  the  interval  0  <  a  <  1),  then 

fg  !  fg  (39) 

Utilizing  these  inequalities,  we  find  the  lower  and  upper  bound  on  the 
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third  moments,  given 


A  and  A 


7  <  7  -  (A  ~  ^ 

1  -  A 


(40) 


We  confine  our  discussion  to  the  lower  bound  since  the  upper  bound  can  be 
treated  by  similar  methods.  To  establish  the  bound,  we  use  the  Cauchy- 
Schwarz  inequality,  Eq.  (36),  with 


f  -  A3/2 


g  -  A1/2 


(41) 


The  realization  of  the  lower  bound  is  obtained  from  Eq.  (37) 

A3/2  .  XA1'2 

whose  roots  are 

A  =  0  and  A  =  X 

- 3 

Therefore,  the  minimum  of  A  is  realized  by  the  probability  density 

P( A)  =  w6(A)  +  (1  -  w)6(A  -  X) 

From  the  given  moments, 


(42) 

(43) 


(44) 


and 


A  =  (1  -  w)X 
/?=(!-  w)X2  =  A  •  X 


Hence, 


7 

X  =  — 


Substituting  Eq.  (47)  into  (45) 

A  A2 

1  -  w  =  y  =  ~ 

X  —j 
fC 

which  completes  the  determination  of  P(A)  .  Using  Eq.  (44), 
P(A)  =  A-  -  6(A)  +  6  { A  -  — 

7  7 


(45) 

(46) 


(47) 


(48) 


(49) 


12 


In  terms  of  the  box  diagram,  we  can  symbolize  this  distribution  as  in 
Figure  8. 

An  entirely  analogous  discussion  of  the  upper  bound  is  obtained  with 

f  =  A  B1/2  ,  g  =  B1/2  (50) 

The  resulting  extremal  distribution  is  of  Figure  8. 

The  statistical  bounds  established  in  Eq.  (40)  are  sufficiently 

~ T 

stringent  to  give  information  on  the  experimental  values  of  B  that  can 

be  obtained  from  Konrad's  paper.  In  Figure  9,  we  show  the  successive 

restrictions  obtained  by  statistical  considerations  when  we  use  the  infor- 

T 

mation  that  B  is  given  and  when  we  use  both  B  and  B  .  Our  estimated 

“T 

experimental  error  is  the  one  discussed  in  the  context  of  B  .  The 

maximum  width  of  1/16  is  the  maximum  of  the  difference  between  upper  and 

~’T 

lower  bounds  in  Eq.  (40).  The  sensiti vity_of  the  bounds  to  errors  in  B 

~~7  ? 

and  B  is  readily  estimated.  For  B  =  BQ  ±0.01,  at  y/(x  -  xQ)  =  0,  is 

6L3  s  ±0.013 

6U3  =  ±0.010 

6W3  2  *0.003  ,  W3  =  0.039 
~7 

If  both  B  and  B  are  varied  as 

B  =  B^  ±0.01 

f?  =  ej  ±0.01 

we  find  comparable  variations  in  the  bounds. 

It  is  of  interest  to  restate  the  bounds  of  Eq.  (40)  in  terms  of  the 
— T 

centered  moment  A'  .  From  the  relationship 

A3  =  A3  +  3A  A*7  +  A-7  (51) 


we  obtain 

<A^<A^  (-1.  ~  (52 

An  immediate  corollary  of  these  bounds  is  the  determination  of  the  allowed 
domains  of  positive  (negative)  skewness.  These  are  shown  in  Figure  10. 
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A  third-order  mixedness,  ,  is  introduced  by  analogy  with  Eq.  (21) 


A  *  B  ’  -  A '  B ' 


U3  -  4 


^3  .  a3-  4 

U3  "  L3 


Clearly,  0  ^  R^  £  1  .  A  plot  is  shown  in  Figures  11a, b  for  two  of 
Konrad's  flows  in  terms  of  mass  fractions.  Rj  appears  to  be  systemati¬ 
cally  asymmetric. 

A  discussion  of  the  fourth  moment  can  be  carried  out  somewhat  analogous¬ 
ly  to  our  previous  discussion  of  the  third  moment. 

The  statistical  bounds  can  be  summarized  in  the  following  results  (the 
bounds  are  proven  in  Section  IV). 

1.  For  given  A  we  have  the  broad  bounds 

A4  <  A*  <  A  (54) 


2.  For  given  A  and  kc  we  have 

~t3 

7>£. 

A 

A4'  <  +  A2  +  (A  -  ?)2.  -  A2 

(  1  -  A)2 
“  T  ~3 

3.  If  the  three  moments.  A,  A  and  A  ,  are  all  given,  then 


~K  T‘ 
IV  >  PT 


ft 3  t  ft2 , 


(A  -  A  A 

~7  -? 

K  -  k 


1?  <  i?  -  l*2  -  fl3)  (5f 

*-7 

A  plot  of  the  successively  shrinking  ranges  of  the  upper  and  lower  bounds 

on  obtained  from  Konrad's  data  is  shown  in  Figure  12. 

The  sensitivity  of  the  stringent  bounds  on  variations  of  the  lower 

~~7  ~ T 

moments  is  approximately  the  same  as  for  B  and  B  . 


6 


"T 

The  relation  between  the  total  moment  B  and  the  Kurtosis  is  obtained 
through  the  Reynolds  decomposition  which  yields 

B4  =  B4  +  6B2  B^  +  4B  jP”  +  b"^"  (59) 

Upon  use  of  (51),  we  find 


Kurtosis  = 


B'4  =  B^~  =  3B4  -  4B  B3  +  6B2  B2 
B^  (?  -  B2)2 


(60) 


A  plot  of  the  Kurtosis  across  the  flow  is  shown  in  Figure  13.  The  Gaussian 
value  of  3  is  nearly  achieved  at  the  center. 
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III.  INCOMPRESSIBLE  FLOW  WITH  VARIABLE  DENSITY 


In  this  section,  we  characterize  the  state  of  an  incompressible  flow 
with  variable  density  by  means  of  appropriate  state  functions.  In  parti¬ 
cular,  the  (mass)  density  and  the  entropy  of  mixing  are  calculated  as 
functions  of  the  mass  fractions,  and  the  relations  among  molar  and  mass 
variables  are  given. 

To  establish  our  notation,  we  recall  that  for  a  single  component  gas 
two  measures  of  density  are  useful;  the  mass  density  p  and  the  molar 
density  ff  .  They  are  related  through  the  number  of  molecules  per  unit 
volume,  n  (number  density) 


p  =  m  n  =  M/T 

where  m  =  mass  of  one  molecule 

M  =  m  x  Arogadro's  number 
=  mass  of  one  mole 

The  gas  constants  are  related  by 

R  =  ®  =  — 
m  M  m 


(61) 

(62) 


(63) 


where  k  is  Boltzmann's  constant.  The  equation  of  state  of  our  ideal  gas 
thus  acquires  two  possible  forms  in  terms  of  molar  and  mass  densities: 


p  =  nkT  =  /TV T  =  pR  T 
m 


(64) 


Correspondingly,  for  a  binary  mixture,  we  have  two  measures  of  concentra¬ 
tion: 

1.  Mass  fractions 


a  = 


P1 

Pi  +  P2 


p2 

Pj  +  P2 


(65) 


with 


and 


a  +  6  ■  1 

0  <  a  <  1  ,  0  <  3  <  1 


(66) 

(67) 
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2.  Molar  fractions 


A '  ^  Tjq 


«2 

B=J^1T2 


with 


and 


A  +  B  =  1 


0  <  A  <  1  ,  0  <  B  <  1 


(68) 


(69) 

(70) 


We  adopt  by  convention  the  correspondences 

1  < — »  a  < - >  A  < - >  heavier 

2  < — B  < - >  B  < - lighter 

The  equation  of  state  for  a  mixture  of  two  ideal  gases  can  be  written, 
assuming  Dalton's  law  of  partial  pressures. 


P  =  P1  +  P2 


=  /f«T  =  (Jfj  +K2)^t 

■  f*T  ■  p(t  +  i^)  ST 


(molar  form) 

(mass  fraction  form) 


(71) 


An  incompressible  flow  with  variable  density  is  characterized  by  the 
assumptions 


p  =  p  =  constant  across  the  flow 
T  =  T  =  constant  across  the  flow 
We  then  see  that  in  molar  variables 

Jf=  J*l+Jiz=  Jf  =  constant  across  the  flow 


(72) 


(73) 


while  in  mass  fraction  variables 


constant 

t)  (v 


t  constant 


(74) 
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We  introduce  the  normalized  density 


where 


_ 2 _  =  1  =  A  a  +  1 

‘'light 

(75) 

„  .V 

plight 

(76) 

** 1 

(77) 

We  can  choose  without  loss  of  generality,  0  <_  A  <  1  .  Thus  for  an 

He  -  mixture,  A  =  6/7  «  0.86  and  for  an  H2  -  0^  mixture,  A  =  15/16 

«  0.94. 

We  observe  that  p*  is  a  rather  steep  function  of  a  for  A  ~  1 
and  thus  p*  is  statistically  independent  from  a  except  for  very  weak 
turbulence.  In  contrast,  p*  is  linear  in  A  and  thus  p*  and  A  are 
equivalent  statistical  variables.  The  relationship  between  the  molar  and 
mass  fractions  is  therefore  nonlinear.  We  have,  in  fact. 


A  =  (1  -  A)  A(a) 


(78) 


B  = 


1  -  a 
1  -  Aa 


(79) 


and  the  inverse  relations 


a 


A 

1  -  A(1  -  A) 


(80) 


^  ‘  (!  '  A)  1  -V-  A) 

A  graph  of  the  relation  A(a)  is  shown  in  Figure  14  to  show  the  increasing 
statistical  independence  between  A  and  a  with  increasing  A  . 

Another  nonlinear  state  function  of  the  mass  fraction  is  the  entropy  of 
mixing.  The  experimentally  observed  pure  Ng  and  pure  He  probability  shows 
that  mixing  is  very  slow  on  the  time  scale  of  turbulence  variations 
(~  integral  scale/root  mean  square  turbulent  velocity).  Further,  the 
probability  profile  of  pure  He  is  approximately  2-3  times  thicker  than  the 
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probability  profile  of  pure  Ng  in  qualitative  agreement  with  the  concept 
that  a  pocket  of  pure  folded  by  a  large  scale  eddy  into  the  He-rich 
region  should  be  mixed  by  molecular  diffusion  about  2.6  times  faster  than  a 
pocket  of  pure  He  folded  into  Ng  =  ^  mN^mHe  =  ^  *  2.65).  We 

show  below  that  choosing  a  probability  distribution  with  low  entropy  of 
mixing  is  in  rather  good  agreement  with  experimental  evidence.  In  this 
section  we  calculate  the  entropy  of  mixing  and  give  some  of  its  properties. 

For  a  single  ideal  gas  we  can  write,  for  the  specific  entropy, 

S  =  Cy  ln(T/T  f)  -  Rm  *n(p/pref)  (82) 

m 

where  Rm  is  given  by  Eq.  (63).  For  a  binary  mixture,  we  have 


S  (mixture)  aSct  +  3S8 


(83) 


The  entropy  density  is  obtained  as  Psspecific-  Substituting  (82)  into 
(83) 

S  .  v+  =  aC.,  S'"  ? - +  &CM  1 - + 

mixt  va  Tref  v3  Tref 

Pf,  Pd 

-  aR  £n  — - 8Ra  Jtn 


a  pref  3  pref 
We  choose  for  the  variable  density,  incompressible  flows  that  we  analyze 


(84) 


Tref  =  Texp  =  constant 

=  5. 

pref  plight(a) 

The  resulting  expression  for  the  specific  entropy  is 


(85) 

(86) 


mixt 


fi 


^n(ap*)  +  jjp  £n(Bp*) 

a  6 


(87) 


which,  using  Eq.  (75),  reduces  to 


M 

S  —  =  - 
mixt  H 


«  «n  a  +  0.  -  -  Si  ♦  1^. 


Aa 


fcn 


1  -  Aa 


(88) 
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The  normalized  entropy  density  is  given  by 


Mr 

s  .  = _ l _  a  in  a  +  a)  &n  U  - a)  -  1  1  Aa  on  ( i  .  Aa\ 

*  mixt  n  1  -  Aa  a  1  -  A  1  -  A  *  u  m) 


The  nonlinear  behavior  of  the  specific  entropy  is  shown  in  Figures  15  and 

16  for  physically  interesting  values  of  A.  In  Figure  17  we  see  how  the 

entropy  density  combines  the  nonlinear  features  of  p*  and  of  S  .  *  . 

*  mixt 

As  a  simple  example  of  the  behavior  of  the  mean  entropy  of  mixing,  we 
observe  that  for  A  =  0  (p*  =  1)  we  have  from  (88) 


-£■  =  -  (A  in  A  +  B  Jin  B)  (9C 

Thus  a  simple  computation  shows  that  the  mean  entropy  associated  with  the 
"nonoverlapping"  probability  of  Figure  1  is  exactly  zero. 


while  the  entropy  of  mixing  for  the  distribution  of  Figure  14  is 


=  -  f-  (A  in  A  +  B  in  B) 
a 


It  is  not  difficult  to  show  that  (91)  is  minimal  for  given  A  while  (92) 
is  maximal  for  given  A  . 

We  shall  use  below  the  distribution  given  by  Figure  7.  In  this  case, 


_  —k  r  in  A.  £n(l  -  A,) 

smix  =  -r  <A-A  >  T^x:+  -XT  3 

a  L  3  3 


T 

Thus  for  the  distribution  of  Figure  7,  Sm^x/(A  -  A  )  is  a  function  of  the 
single  parameter  A^  and  can  be  minimized  analytically. 


IV.  GENERAL  METHOD  FOR  THE  DETERMINATION  OF  STATISTICAL  BOUNDS 


In  this  section,  we  first  outline  a  general  technique  to  obtain 
statistical  bounds  that  are  sharp  (i.e.,  "best  bounds")  and  to  establish 
the  domain  of  moments  for  which  they  are  realized.  We  then  give  several 
families  of  useful  bounds. 

The  essence  of  the  method  is  the  close  relation  between  the  eaualitv 
condition  in  the  Cauchy-Schwarz  theorem  and  the  characteristic  equation 
for  a  discrete  probability  distribution. 

Consider  a  discrete  probability  distribution  of  general  form 

n 

P(A)  =  ^  ei6(A  -  A.)  (94) 

i=l 

Using  Dirac's  identity 


(A  -  Ai)6(A  -  A.)  =  0 


(95) 


we  see  that 


TT  (A  -  A.)P(A)  -  0  (96) 

i=l 

so  that  the  algebraic  equation 

n 

♦  (A)  -  TT  (A  -  V  =  0  (97) 

i=l 

determines  the  locations  of  the  masses  of  P(A)  .  We  call  <P  =  0  the 
characteristic  equation  for  (94).  We  can  actually  determine  the  A^'s  as 
functions  of  the  moments  as  follows.  Multiply  (96)  by  any  function  F(A) 
and  integrate  from  0  to  1,  obtaining 


(f(a)TT(a  -  Ai^  =  0 


(98) 


If  we  choose,  in  particular,  F  =  A°=l,  A,  A  ,  ....  A"’1  ,  we  obtain  n 
linear  equations  for  the  invariants  of  the  matrix 


Lij  '  A(i)6ij  (no  sum  on  i  ) 


(99) 
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For  example,  for  n  =  3  ,  the  linear  equations  are 


where 


A3  -  Tlj  +  AI2  -  I3  =  0 

~T  “7  ~7 

Iv  -  +  A^I2  -  AI3  =  0 

tt  -r  "T  ~7 

A3  -  A^Ij  +  AjI2  -  A^I3  =  0 


A1  +  a2  +  a3 

AlA2  +  AiA3  +  A2A3 
A1A2A3 


( 

(101 

(102 


(104; 


(105) 


By  writing  Eqs.  ( 100 ) - ( 102 )  as  a  matrix  equation 


A2  -  A  1 

h 

’  ?  ’ 

7-7  A 

h 

= 

7 

7-7  7 

h 

7 

- 

_  _ 

_  _ 

(106) 


we  can  obtain  the  Ik  's  in  terms  of  the  Cramer  determinants  of  (106). 


Thus,  with  obvious  notation 


N, 


I  =  — 
k  Nq 


(107' 


It  follows  then  quite  generally  from  Eq.  (97)  that  the  A^  's  are  the 
n  roots  of  an  n'th  order  polynomial  whose  coefficients  follow  from  (107) 


N  A" 
o 


»n-L 


(-I)V 


=  o 


(108 


The  locations  A.,  of  the  masses  are  then  determined  by  establishing  that 
N..  >  0  and  Nn  c  NQ  .  NQ  is  always  positive  by  Hankel's  inequality. 
This  is  established  as  follows: 


2 


Substitution  of  Eqs.  (112)  and  (113)  into  (111)  yields 

A2"  >  l(pF,  A1,  ...,  A2n_1  ) 

To  obtain  the  lower  bound  on  A cu  ,  the  upper  bound  on  PC  ,  and  the 

upper  bound  on  A2 0+2  ,  we  multiply  (112)  and  (113)  successively  by  A"2  , 

B2,  (AB)2  .  For  the  lower  bound  on  A^n  ,  we  return  to  the  forms  (112) 
and  (113). 

As  a  simple  example  of  the  method,  we  establish  the  least  upper  bound 
~T  -  2  3 

on  A  for  given  A,  A  ,  and  A  .  Since  we  require  the  upper  bound  on 
an  even  moment,  we  consider  the  distribution  of  Figure  18.  The  character¬ 
istic  equation  is 

A(A  -  1 ) (A  -  A3)  =  0  (114) 

averaging 

(115) 


<^A(A  -  1)a)>  -  A3<(A(A  -  1))>  =  0 
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Thus, 


7 


To  show  that  <  1  ,  we  observe  that  <^A(A^  +  1  -  2A)^>  j>  0  .  The 
location  of  the  "middle"  mass  is  thus  determined  without  knowledge  of  the 
weights.  To  find  e  we  note  that 


Thus, 


A(A  -  1)  =  eA3(A3  -  1) 


~7  -  3 

(A2  -  A) 


(7  -7)  {7  -  2?  +  A) 


(117) 


(118) 


the  proof  that  e  <  1  is  analogous  to  the  one  for  A3  <  1  .  Having  estab 
lished  the  realizability  of  the  limit,  we  construct  the  Cauchy- Schwarz 
inequality  with 


f  =  A^l  -  A)^A 


(119) 


9 


k2{  1 


(120) 


and  find  after  simplifications  that 


A 


4 


(121) 


We  note  that  the  proof  given  establishes  Eq.  (121)  as  the  "best  bound"  for 

“  l 

tC  in  the  following  sense.  The  right-hand  side  of  (121)  is  the  least 
upper  bound  on  7  for  all  probability  distributions  of  given  A,  7  , 

and  7  . 

An  important  shortcut  for  the  odd  moments  is  obtained  from  the 
"renormalization  theorem."  Given  that  we  have  established 


f2(A)D1(A)  da 


1  ^ 

g2(A)D1(A)  dA  > 


(122) 
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it  follows  that 


(123) 


The  proof  is  simple.  Since (122)  holds  for  all  ,  we  can  choose  D1  of 
the  special  form  for  any  D2  .  (The  passage  from  (123)  to  (122)  is 

not  legitimate,  in  general.) 

We  now  consider  several  important  families  of  inequalities.  First, 
consider  a  variable  a  ,  0  <_ a  <  1  ,  and  its  conjugate,  0  =  1  -  a  , 

i.e.,  either  a  mass  or  a  mole  fraction.  Applying  Cauchy-Schwarz  to 


f  =  a  ,  g  =  1 


(124) 


we  obtain 


1  7 


T 

a 


(125) 


where 


1 

D(a)  da  . 


If  D  is  a  probability  density,  1  =  1  .  The 


presence  of  the  1  allows  application  of  the  renormalization  procedure. 
The  Cauchy-Schwarz  limit  for  (125)  is  a  =  A  ,  i.e.,  a  single  mixed  cell 


P(a)  =  <5(a  -  a) 

Renormalizing  (125)  by  a  we  find 


a 


The  Cauchy-Schwarz  limit  is  a^2  =  Aa1^2  ,  i.e.,  a  =  0,  a  =  A 
Therefore,  at  equality, 

p(o)  =  £&(a- £)  + 

TV  a  ~7 

a  '  '  a 

Renormalizing  Eq.  (125)  by  0  ,  we  find 

_ 2  ~T  - 

alT  <  a  0  *  B 


(126) 


(127) 

T 

a 

a 


(128) 


(129) 
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or,  eliminating  6  , 


~T  _  ~7  (a  -  a  ) 

a  <_  a  -  - L- 

1  -  a 

The  Cauchy-Schwarz  limit  for  Eq.  (127)  is 


(130) 


with  roots 


aB  2  =  XB 2 


6  =  0,  a  =  X  = 

6 


(131) 

(132) 


i.e.,  a  pure  a  and  a  mixed  cell.  Therefore,  at  equality 


T 


P(a)  =  °?l.  6  (a  -  1 


T  -2 


— - 
1  -  2a  +  a 


1  -  a  /  -  9-  2 

'  1  -  2a  +  a 


6(a) 


(133) 


Renormalizing  (125)  by  a6  ,  we  find 


TT 


“T“  _ 


a  6  <  a  6  •  aB 


(134 


or,  eliminating  B  , 


T  “X 


~J  „  ~5  (a  -  a) 

a  <  a  -  - L 


~1 
a  -  a 

The  Cauchy-Schwarz  limit_was  discussed  following  Eq.  (114). 

4 

The  lower  bound  on  a  requires  two  mixed  cells.  We  find 

T .  “?2  .  (?-  ZJ)2 

a  >  a  +  - L— 


T  -2 

a  -  a 


(135 


(136) 


The  a  inequalities  and  some  of  the  higher  ones  are  summarized  in  the 
trees.  Figures  18  and  19.  For  all  cases,  the  equality  corresponds  to  a 
unique  discrete  distribution.  This  situation  persists  when  we  include 
"skips."  These  are  most  readily  derived  from  a  generalization  of  Eq.  (125) 
which  is  obtained  from  Holder's  inequality 


I"-2 


n-1  -n-1 

a  >  a 
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(137) 


We  observe  that  for  n  =  3,  Eq.  (137)  becomes  (125).  We  readily  find, 

3 

following  the  tree  of  Figure  20,  the  bounds  on  a  given  only  a  ,  i.e  , 


-3,1- 
a  <  a  <  a 


(138) 


The  equal  sTgn  of  the  lower  bound  is  realized  by  (126)  just  like  the  lower 

2  4  _ 

bound  on  a  .  Similarly,  the  bounds  on  a  given  a  and  a  are 


•  ,1,1 

~2  —  a  —  a 


[a  -  ?) 
(1  -  a)4 


(139) 


and  they  are  realized,  respectively,  by  (128)  and  (133).  In  Appendix  3 
we  prove  a  strong  version  of  this  result. 

The  inequalities  for  p*a  are  most  readily  obtained  by  observing  that 
renormal izing  Eq.  (125)  oy  the  (positive)  quantity  p*  gives  us 


P*  '  P*a  1  P*a ' 


(140) 


Using  the  equation  of  state  and  the  relations  given  in  Appendix  4,  we  can 
recast  (140)  either  as  a  lower  bound  on  p^  (given  a) 


1  -  Aa 

or  as  a  lower  bound  on  p*a  (given  a) 


(141) 


P*a  > 


A  L  1  -  Aa 


(142) 


Following  the  tree  of  Figure  21  or  the  (simpler  but  less  obvious)  alterna¬ 
tive  tree  of  Figure  22  and  using  the  equation  of  state,  we  readily  obtain 

_  _  - 2" 

the  lower  and  upper  bounds  on  p*a  for  given  a  and  a  as  well  as  the 

_  _  3  _ 

upper  bound  on  p*a  for  given  a  a  and  a  .  The  lower  bound  on  p*a 
-  T  3 

for  given  a,  a  and  a  initiates  a  new  tree  and  is  relatively  difficult. 
We  thus  give  it  in  detail.  We  show,  in  fact,  that  given 


T  -  -1 
1  ,  a  »  a 


we  have  the  upper  bound  on  a' 


(143) 


(144) 
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where 


~r 
~J  e*a 

rt  —  -  ■  — — 

[?(1  -  A2?)  -  a(a  -  A?)  ] 

-  ~2  T2  ~3 

-  a  a  -  Aa  +  a° 

u*  - 

A 

p*a  (1  -  Aa)  -  a 

7 

Clearly  the  (least)  upper  bound  on  a  is  equivalent  to  the  (greatest) 

_  p  o' 

lower  bound  on  p*a  (given  a, a  and  a  ). 

We  choose 


To  obtain  the  Cauchy-Schwarz  inequality,  we  square  f  and  g  ,  obtain  the 
product  fg  ,  and  average.  The  result  is 


+  2c 
+  c2 


rranged  as  coefficients  of  cn  to  facilitate 

n  of  state  relations  of  Appendix  3,  we  obtain  from 
ancellations ,  an  inequality  which  is  linear  both  in 
slight  variant  in  the  choice  of  f  and  g  was 
his  goal ) : 


f  =  P*  (a  +  c) 


9  =  P*  a 


(145) 


(146) 


p*a  -  p*a  •  p*ct 


(1  +  Ap*a)p*a  -  p*a" 


>  0 


(147) 


"  P*a 

-  T  "Z2  -  T~l 

,  a  a  a  a  a 

.?■> 

~z — 2 - r 

A  hr  hr 

(f  -  a),  f]  >  0 


(148) 


To  ensure  that  the  inequality  yields  the  least  upper  bound  on  (over  the 

_ 

probability  densities  of  given  1  ,  a  ,  a  ,  and  p*a  ),  we  rewrite  (148) 


(149) 


F(c 


)  > 


and  minimize  F  with  respect  to  c  .  We  have 


F(c)  » 


p^ci  T  T  -  T 

*  a  a  a  a 

~7T  A  TT“  .3 

A  A  A 


+  2c 


~ 


a  .  ~ \  .  a 
-  t  +  a  I  +  —k 

A 


c2[p*a(i-  a)  -  f]  (150) 


The  curvature  of  F  is  positive  by  virtue  of  (142)  and  the  first  deriva¬ 
tive  vanishes  at  c  =  c 


_  _2 

,  p*a  (a  -  Aa  )  -  a 

s =  — : — r 

p*a  (1  -  Aa)  -  a 


(151) 


Substituting  this  value  of  c  into  (149)  we  obtain,  after  some  algebra, 
3  — T 

a  <■  a*  as  given  in  (144). 

We  now  examine  the  equal  sign  condition,  i.e.. 


f  =  Ag  (152) 

Using  (145)  and  (146),  we  have 

p  *  (a2  +  c)  =  Ap*  a  (153) 

Since  p*  has  no  zero  for  a  between  0  and  1,  Eq.  (153)  gives 

a2  -  aA  +  c  =  0  (154) 


We  now  determine  A  and  c  by  averaging  the  characteristic  equation 
(154)  and  also  multiplying  by  a  and  averaging: 


7 

~T 

a 


aA  +  c  a  0 
~7 

aX  +  ac  =  0 


(155) 

(156) 


We  have  as  solutions  of  (155)  and  (156): 


(157) 
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1 


c 


-  "7 


_  a  a 


T 

g 


(158) 


It  is  clear  from  (154)  that  the  probability  density  that  realized  the 
equality  is  the  two-mixed  cells  distribution: 


P(g)  =  e^6(g  -  g^)  +  £2(5(a  "  “2) 

where 


al  +  a2  =  ^ 


(159) 


(160) 


ala2  =  c 


(161) 


as  given  by  Eqs.  (157)  and  (158).  In  order  to  complete  the  realizability 
proof  of  the  equal  sign,  we  must  show  that  (159)  has  all  four  parameters, 
,  £2  .  aj  and  between  0  and  1. 

We  first  discuss  the  locations  and  prove  that 


al  +  a2  —  0 

(162) 

ala2  >  ® 

(163) 

and 

ala2  £  1 

(164) 

Clearly,  Eqs.  (162),  (163)  and  (164)  imply  that  g^  and  g^  are  between 
zero  and  one.  To  prove  these  inequalities,  we  multiply  the  characteristic 
equation  for  (159)  by  an  arbitrary  function  F  assumed  similar  to  a  .  We 
have 


^T(g^  +  g£)  +  Fa  ^2  =  0 


and  (F  =  1) 

a  -  aUj  +  a2)  +  =  0 

Equations  (165)  and  (166)  determine  g^  +  g2  and  g^2 


(165) 


(166) 
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(167) 


al  +  a2  " 


F  aL  •  j  -  F  Z 
Ta  •  1  -  F  •  a 


and 


-  — T  —  T 
_  a  F  a  -  Fa  a 
a  a  -  — - = - 

1  c  Fa  •  1  -  F  •  a 


By  using  Tchebycheff 1 s  inequality,  we  find 


Fa  •  1  >  F  a 


because  F  is  similar  to  a  .  Also 


— 2  -  ~  ~2 
Fa  •  1  >  F  a 


because  F  is  similar  to  a  and  renormalizing  (169)  by  a 


— 2  -  —  T 

Fa  •  a  >.  Fa  *  a 

This  completes  the  proof  of  Eqs.  (162)  and  (163). 
To  prove  that  a^  <  1  ,  we  need 

.  — 2  —  T  —  -  ~  - 

a  Fa  -  Fa  a  <  Fa  •  1  -  F  •  a 


(168) 


(169) 


(170) 


(171) 


(172) 


We  rewrite  (172)  as 


<\F  ip(a)^>  >_  0 


with  \J>  given  by 


~2  2  - 

\p(a)  =  a  •  1  -  a  +aa  -  a  a 


But  ip  is  similar  to  a  because 


i|)'  (a)  =  I  +  a  -  2a  a 


is  linear  in  a  and 


ip'(0)  =  I  +  <7  >  0 


(173) 


(174) 


(175) 


(176) 
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so  that 


We  also  note  that 


V(l)  =  (1  -  a)  +  a'"  >  0 


(a)  >0  for  0  <_  a  <_  1 


ip  =  0 


Applying  Tchebycheff  to  (173) 


Fif/  >_  F  tp  =  0 


q.e.d. 


We  have  completed  the  proof  that  and  are  between  0  and  1. 

We  note  that  a,  +  a~  and  a,a0  can  be  expressed  either  in  terms  of 

— g-  1  _ £  1  t 

a  or  in  terms  of  p*a  .  This  shows  that  the  two-mixed  cells  distribu- 

—  -  -  "7  “3  -  _  ~T  _ 

tion  realizes  either  (1 ,  a,  a  ,  a  )  or  (1 ,  a,  a  ,  p*a). 

We  now  prove  that  the  weights  Cj  and  Ed-  (159)  are  between 

zero  and  one.  By  introducing 

=  a  -  6 ^  (18 


we  find  that 


imply 


a2  -  a  +  &2 


ct  =  0 
T  =  1 


E1  "  «1  *  «2  ’  '  4l  +  «2 

A  simple  calculation  further  shows  that 


T.  .  , 

a  o  jo  2 


771  4i 


Since  a'  is  given  (and  positive),  6^  and  62  have  the  same  sign.  I 
follows  then  from  (185)  that  Cj  >  0  and  >  0  •  Since  their  sum  is  1 
both  e,  and  e„  are  between  zero  and  one. 


2 

Another  important  family  of  inequalities  is  the  one  pertaining  to  p*  . 
Clearly,  from  Cauchy-Schwarz 

1  •  >  pT2  (187) 

~~2  _ 

This  is  a  lower  bound  on  p*  for  given  p*a  .  The  tree  of  Figure  23  yields 
the  relevant  bounds  except  for  the  lower  bound  on  pf  for  given  a,  a  , 
p*a  .  We  proceed  to  prove  this  bound  separately. 

We  choose 

g  =  -p*^  (188) 


f  =  <t>  +  9  =  P*(ot  -  »2)  "  P*^  (189) 

where 

4>  =  (P*  "  P*1)(a  -  ot2)  (190) 


P*1  =  1  -  Aa} 


(191) 


Thus,  Eq.  (190)  is  the  characteristic  equation  for  two  mixed  cells  with  one 
of  the  singularities  expressed  in  terms  of  p  .  We  use  the  equation  of  state 
relations  of  Appendix  3  to  eliminate  products  of  p*  and  a  (to  a  power). 
After  some  algebra,  we  find 


f2  2 
f  “  P* 


(° 2  -  i)  +  p*f  (pn  +  1>(“z  ■  i)* 


,  2  2 
+  a  p* 


1  +  I  apl*  + 


2p 


*la2 


(192) 


g  =  (p* 


la2) 


(193) 


(194) 


By  using  the  Cauchy-Schwarz  inequality,  we  find 


"7/  l 

\2  —2 

f  1\ 

A  7  2 

-  2  . 

^  1 

p*(a2'A 

)  +  p*  A  (p*l  +  1) 

(a2  ”  A  / 

+  a  p*j 

+  a  A  P*1  +( 

u2  a2  p  yj 

i  [p*(}‘“2)'  “p*l  '  S. 


(195) 


The  fact  that  the  inequality  represents  a  lower  bound  on  p*  follows  from 
the  fact  that  its  coefficient  on  the  large  side  is  positive.  The  fact  that 
the  bound  is  the  greatest  lower  bound  for  given  ci  ,  c?  ,  and  p*a  follows 
from  the  realizability  of  the  two-mixed  cells  distribution  which  was 
obtained  above,  Eqs.  (152)  through  (186). 

The  last  family  of  inequalities  that  we  consider  are  bounds  on  the  mass 
fraction,  given  the  mole  fraction.  From  the  tree  of  Figure  24,  several 
follow. 
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V.  MODEL  SELECTION 


In  this  section,  we  first  describe  our  process  for  selecting  a  model 
probability  density.  Having  chosen  a  model  with  a  free  parameter  that  is 
adjusted  to  minimize  the  average  entropy  (density),  we  compare  several 
features  of  the  model  against  the  data  reduced  and  presented  by  Konrad. 

In  selecting  the  model,  we  find 

1)  At  the  edges  of  the  mixing  region,  the  experimental  probability 
distribution  is  sufficiently  skewed  to  rule  out  symmetric  models. 

2)  The  experimental  probability  for  finding  pure  N2  and  pure  He  is 
sufficiently  large  to  rule  out  models  that  do  not  have  strong  spikes  at 
molar  fractions  close  to  one  or  zero. 

3)  We  find  that  the  third  moment  is  rather  insensitive  to  the  selection 
of  a  model  providing  the  model  is  statistically  correct. 

We  conclude  that  the  most  satisfactory  model  distribution  is  one  for 

which  the  average  entropy  density  is  minimized. 

With  regard  to  the  skewness  factor,  we  consider  first  a  square  hat  in 

A  (Figure  25).  The  probability  density  has  three  parameters  and  it  is 

—  T 

thus  fully  specified  by  T  ,  A  ,  A  .  The  moment  equations  are 

1 

P(A)  dA  =  H(R  -  L)  (196) 


A-  V 

(197) 

A2  =  -j  (R2  +  RL  +  L2) 

(198) 

with  solutions  for  the  edge  locations 

R  =  A  + 

(199) 

L  =  A  - 

(200) 

where 

(201) 

and  finally  for  the  height 

(202) 

2/  3  A'2 
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The  fit  to  Konrad's  equal  mass  flow  is  shown  in  Figure  26.  Unacceptably 

large  probabilities  are  found  for  negative  molar  fractions.  The  origin  of 

the  difficulty  is  the  large  skewness  of  the  edge-flow  distributions.  The 

point  is  illustrated  in  Figure  27  by  a  simple  triangle  distribution  that 

cannot  be  fitted  by  a  legitimate  square  hat.  The  region  of  moments, 

—  jf 

A,  A'  ,  for  which  a  square  hat  is  statistically  acceptable  is  readily 
obtained  from  R  <  1  and  L  >  0  .  This  region  is  shown  in  Figure  28, 
together  with  the  region  allowed  by  the  less  restrictive  but  still  symmetric 
model 

A'3  =  0  (203) 

The  various  boundary  curves  are  given  by  the  following  functions  of  A  : 

1  :  A2/3 

2  :  B2/3 

3  :  2 A2  (204) 

4  :  2B2 

5  :  A(1  -  A)  =  AB 


In  Figure  29,  we  show  the  value  of  the  ratio 


where 


r  = 


^model _ 

^experiment 


Q  = 


(205) 


(206) 


The  value  of  Qy  corresponds  to  the  upper  bound  for  Q  given  A  and 

A'  ;  similarly  Q.  corresponds  to  the  lower  bound  for  Q  ,  given  A  and 
"  2  *“ 

A  .  The  graph  emphasizes  the  fact  that  towards  the  edges  of  the  mixing 
region  any  symmetric  model  would  be  in  disagreement  with  experiment. 

A  model  with  much  less  symmetry  and  some  promise  of  describing  a  flow 
with  low  mixing  is  defined  by 


0  =  (AB)'  =  AB 1  +  A'B  +  A'B'  -  A 1 B ' 


(207) 


Equation  (207)  is  a  quadratic  for  A'  once  B’  is  eliminated  (B‘  =  -A') 
and  the  resulting  probability  density  is  given  by 
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where 


and 


P(A)  =  e^A  -  Ax)  +  e26(A  -  Ag ) 


Al  = 


1  -  D  „  _  1  +  D 

2  ’  m2  2 


D  =  A  -  4AB 


We  note  that  for  this  model 


(208) 

(209) 

(210) 
(211) 

(212) 


A 


2 


(213) 


Furthermore,  Eq.  (207)  is  exactly  equivalent  to  the  requirement 
A2B2  =  min  (for  fixed  A,  /?  ) 


(214) 


In  Figure  30,  the  values  of  e1  and  of  e2  are  superimposed  on  Konrad 
probability  densities  and  are  seen  to  follow  the  patterns  reasonably  well. 
The  model  is  not  acceptable  in  view  of  the  experimental  values  for  pure  He 
(Fig.  31)  and  pure  N2  (Fig.  32).  The  model  given  by  (208)  yields  zero 
values  for  both  probabilities. 

In  view  of  the  experimental  probabilities  for  pure  species,  the  model 
given  by  in  Figure  33  offers  a  chance  with  both  pure  cells  built  in. 
Comparison  with  the  probability  density  is  shown  in  Figure  32  and  is 

tolerable.  Comparison  with  the  pure  species  data.  Figures  31  and  32,  shows 

however  that  the  weights  of  the  pure  cells  are  too  weak. 

Since  the  entropy  of  mixing  is  clearly  quite  low  in  the  flow,  the  model 

of  Figure  7  was  made  determinate  by  minimizing  the  entropy  (density)  of 
mixing  with  respect  to  A^  .  The  resulting  model  is  defined  as  follows. 

The  weights  are  those  given  by  Eq.  (30)  and  the  bounds  on  A^  are  those  of 
Eq.  (31).  For  this  model,  the  quantity 


n 


1-7 


(215) 


is  sketched  in  Figure  34.  The  quantity  n  is  even  in  A3  with  a  minimum 
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of  approximately  2.76  at  A3  =  0.5  .  A3  is  therefore  determined  as 

Then 


follows.  We  call  the  bounds  of  A3  L  and  U 


T 

A  -  A* 


U  =  4- 


1  -  A 
T 
A 


Minimization  occurs  as  follows: 


1 


If  U  <  j-  ,  then  A3  =  U 


If  L  <  j  <  U  ,  then  A3  =  ^ 


1 


If  ^  <  L  ,  then  A3  =  L 


(216) 

(217) 


(218) 

(219) 

(220) 


The  resulting  weights  are  plotted  in  Figure  35  and  show  a  good  pattern 
relative  to  the  probability  densities;  furthermore,  the  pure  species  plots 
of  Figures  31  and  32  are  in  reasonable  agreement  with  experiment.  The 
fact  that  the  model  distribution  is  somewhat  He  rich  is  in  qualitative 
agreement  with  the  low  mixing  allowed  by  large-scale  transport  and  with  the 
data  on  pure  species. 

By  way  of  contrast  on  this  point,  in  Figure  36  we  show  the  entropy 
maximized  for  the  same  cell  structure.  It  is  in  poor  agreement  with  the 
pure  species  probabilities. 

In  Figure  37,  we  give  the  r  ratio  of  Eq.  (205)  for  the  minimum  entropy 
model.  Although  the  agreement  with  experiment  on  this  point  is  low  (~50%) , 
the  opposite  extreme  of  maximal  entropy  is  roughly  as  good  on  this  count. 
Statistical  consistency  seems  to  be  all  that  is  needed  to  obtain  a  reason¬ 
able  r  .  This  point  will  be  taken  up  again  later  in  connection  with  the 
flows  with  large  mass  differences. 

In  Figures  38,  39  and  40  we  give  the  average  entropy  density,  its 
fluctuations  and  its  second  moment  to  give  an  idea  of  the  actual  flow 
values  and  their  model  values. 

We  now  turn  to  the  discussion  of  the  model  chosen  for  large  mass 
differences.  In  this  case  we  have  at  our  disposal  both  the  mass  fraction 
and  the  density  moments  (equivalently  the  mass  and  the  molar  fractions). 

We  have  thus  chosen  a  four-cell  distribution  with  one  free  parameter. 


_  2 

Figure  41.  We  give  T  =  1,  a  ,  a  ,  p*a  and  p*  .  We  note  that  for 
-J  _  -o 

given  a  ,  a  and  p*a  ,  P*m^n  requi res  two  mixed  cells  (a  discrete 

-jr 

distribution)  and  similarly,  P*max  requi res  two  pure  and  one  mixed  cell 
(again  a  discrete  distribution).  The  algebraic  solution  of  the  moment 
equations  can  be  summarized  as  follows.  The  weights  are  given  in  terms  of 
parameters  and  u*  as  follows: 


e2 

e3 

e4 


=  r 

u  -  u* 

ci 

1  -  u 

-  r 

ui;  -  u 

C2 

u 

=  r 

u3  ‘  u 

C3 

U  -  CC 

=  r 

“4 

L4 

u(l  -  u 

(221) 

(222) 

(223) 

(224) 


The  singularities  of  the  denominators  correspond  to  one  of  the  mixed  cells 
degenerating  with  a  pure  one  or  the  second  mixed  cell. 

The  parameters  and  uT  are  in  turn  given  by 

,  .a  -  p^a  (1  -  AaJ 

C  =  - 1-  (225) 


a  -  1  -  [l  -  p*a  (1  -  A)  (1  -  Aa3) 
Aa-, 


(226) 


(1  -  Aa.)  p*a  (1  -  A)  -  a 

r  =  - - =L  (227) 

Aa3(i  -  a3) 

C4  =  5  (  }  -  1  -  a3  )  -  (1  -  Aa3)  +  7  (228) 


and 


1  -  Aoo  _ 

a - — —  (p*a  -  a) 

a  -  (1  -  Aag)  P*a 


(229) 
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(230) 


1  -  Act 


2a  -  1  -  a  - 


~  [  a  -  P*a  (1  -  A)] 


1  -  a  -  (1  -  Aa3)(l  -  p*a  [l  -  A]  ) 


~1  1  -  A  , -  -x 

a - —  (p*a  -  a) 

a  -  (1  -  A)  p^a 


(231) 


4  A 


(232) 


The  concentrations  are  related  by  the  following  formula  which  is  derived 
using  two  s  cells  and  two  a  cells  in  the  characteristic  equation 

(a2(p* " p*  [l  +  lh]  +  T^a))  ‘  a3^(p*  •  p*[  1  +  TTa  ]  +  lh)) 

(a(p *  -  p*  f1  +  tVJ  +  rh))  ~  a3^*  - p*  f1  +  ih  j  +  ih) 

(233) 


The  minimim  of  the  mean  entropy  occurs  when  =  0  leaving  a  "helium 
rich"  model .  If  we  set 

u  =  u*  (234) 

using  Eqs.  (233)  and  (229),  we  obtain  a  quadratic  equation  for  a^  .  Its 
two  roots  correspond  to  the  presence  of  two  mixed  cells  and  the  interchange- 
ability  of  cell  3  with  cell  4.  The  model  is  thus  fully  determinate. 

We  now  compare  the  model  with  the  two  flows  of  i  /  0  given  by  Konrad. 
We  give 

1)  the  ratio  T  of  Eq.  (205)  (Figs.  42  and  43), 

2)  the  probabilities  of  pure  species  (Figs.  44  and  45), 

3)  the  locations  and  strengths  of  the  Dirac  functions  superimposed 
on  the  Konrad  densities  (Figs.  46  and  47). 

On  all  three  grounds,  the  chosen  model  fares  quite  well.  By  contrast,  we 
show  in  Figure  48  a  maximal  entropy  model.  As  was  the  case  for  A  =  0  , 
the  resulting  model  is  substantially  worse  than  the  minimum  entropy  model. 

As  a  final  indication  of  the  stringency  of  our  statistical  bounds,  we 
give  in  Figures  49  and  48  the  bounds  on  the  (mass  fraction)  moment  ratio 
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with  Q  =  ag  .  The  bounds  are  of  the  same  order  of  magnitude  as  those 

given  in  the  molar  plots.  The  successive  tightening  of  the  bounds,  given 

~ T 

2,  3  and  4  moments,  is  clearly  indicated.  The  bounds  on  a  that  were 
used  in  the  plot  were  not  fully  included  in  the  theoretical  discussion. 
They  are  recorded  below  for  completeness. 


L  <  o?  <  U 


(235) 


~7  _ 

For  given  a  ,  a  ,  p*a 


l.7- 


(a  -  7)  [  |  -  p*a  1  l  A  -  a  +  7_ 


a  -  p*a  (1  -  A) 


-  r 

P*a 


7(1  -  A2  a2  -  a(a  -  A?)j  -  a  a2  -  Aa2  +  ci3 


U  = 


A[  p*a  ( 1  -  Aa)  -  a  ] 


(237) 


(237) 
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APPENDIX  1.  CONVEXITY  OF  MOMENT  DOMAINS 

Let  rrij  ,  in2  e  Mn  be  given. 

For  X  €  (0,1)  we  show  that 

in  =  Xmx  +  (1  -  x)in2  €  Mn  (1.1) 

Write  and  m2  as 

~  (Cq  »  c j  ,  . . .  ,  (1-2) 

~  (^o  *  di  ’  *  ^n-l^  (1-3) 

Then  by  computation  (i.e.,  substitute  (1.2)  and  (1.3)  into  (1.1)) 


m  -  XCg  +  ( 1  -  X )  dg  ,  Xc^  +  (l*X)dj  ,  ...  ,  Xc^_2  +  (1-X)d^_^ 

Since  ir^  €  Mn  ,  there  is  a  non-negative  measure,  Pj  ,  whose 
moments  are  the  components  of  : 


-► 


where 


(Cq  »  Cj  ,  •*.*Cp_|)  "(I,  CL  ,...,  (X  ) 


In  C  np 
a  =  |  aPl 

J  n 


(a)  da 


Similarly  there  is  a  P2(a)  such  that 

m  -  /  A  A  A  \  -  /2  7  2-  2  n-1. 

^2  ~  (do  *  »•  •  •  »dn-r  '  ^  *  a,.. .a  ) 


Consider  a  new  non-negative  measure: 


P  =  XP:  +  (1  -  X)P2 
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Then, 


(1 »  a,  •  •  . ,  a11  *)  = 

=  (Vi  +  (1  -  A)2!  ,  A^a  +  (1-A)2ci  .....  A1^-1  +  (1-A)2an_1) 

=  (acq  +  (1  -  A)d0  ,  AC1  +  (1  -  A)dj  .....  Acn_1  +  (1  -  A)dn_1))=  m 

Therefore,  the  point  m  is  realized  precisely  by  the  measure  P; 
that  is,  iri  €  Mn  . 


PWP 


-2 


a  , 

g  given,  sg  =  - - == 

g  -  Ag^ 

8 

S  =  1  +  ASg 

T 

g 

g 

-y  T  -2 
„2  _  g  -  g 

=  min 


~J  -2 

g  g 

T  ~2 

g  g 


-4 

g 


t  t  [-  rf 1 

a  g  [g  -  Aa  J 


_  "T  . 

6.  g  ,  g  given  ,  Auxiliary  3  =  1  -  g 


g 

a8 

3 

7* 

82 

7 

7 

7? 


T 


Sg  = 


3  =  1  -  2g  +  g' 


H3=g-7 


gS  -  A(g  -  7 ) 


(1  -  A)^l  -  ci  -  A  [a  -  <?  Jj 


=  max 


s  =  1  +  ASg 


7  -  7  -  f 


x2 


-M  +  ^^-i 


3 


i  -  AJ  T 
/  3 


1  -  A2£ 

e  j 


“7  _  -y 

7.  g  ,  g  ,  sg  given,  s  is  maximal 


g 

8 

°3 

C1 

e2 

e3 

(a)  «3  - 


sJ-S 


g  -  (1  -  A)sg 


2 

g  -  ga. 


T 


.  T 

3  ,-g-g  _g-g 

el  "  1  -  g3  *  e2  "  1  '  a  "  g3  *  e3  "  a3(l  -  a3) 


(c)  1=1+  ASg 


(d)  7 


f-ei 

i  =  l  (1-Ag.)‘ 
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APPENDIX  3.  RELATIONS  IMPLIED  BY  THE  EQUATION  OF  STATE 
(REDUCTION  FORMULAE) 

The  following  relations  are  given  as  instantaneous;  they  also  hold  as 
averages. 

1 .  p*  -  Ap*a  =  1 


2. 


3. 


4. 


2  _  1  1  _  1  1  1 
p*a  A  p*a  "  A  a  '  a2  p*  "  A  a  “  a2 


3  _  1  1  1  2  _  P*  a  a  1 

P*a  -  p  P*a  ■  2  a  ~  A  a  ~  3  2  ~  A  ~  3 

A4  hr  h  A  h 


4  1  a  1  2  1  3 

P*a  =  -T  P*«  -  -T  “  “o’  “  -  a  a 

A-3  A  hr 


1  1  12  13  1 

p^  a  a  -  _a  -t- 

AAA  A 


2  12  1 
P*a  ~  A  P*  "  A  P* 


P2.a2  (p.  -  l)2 

hr 

P*ct3  =  4  (P*  -  l)2  -  \  (p*  -  1)  + 

A'3  A4  A4 

P2a4  “  “V  (P*  "  l)2  -  (p*  '  ^  +la+T 

hT  A  h  hr 
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APPENDIX  4.  APPLICATION  OF  INSTANTANEOUS  BOUNDS  TO  ENSURE  CORRECT 
STATISTICS  FOR  SOLUTIONS  OF  THE  RATE  EQUATIONS 


We  consider  a  simple  example  of  second-order  chemistry  in  which  only 
two  species  participate  (although  the  reaction  collision  requires  three 
bodies  as  dictated  by  energy  and  momentum  balance).  The  reaction  chosen  is 


0  +  0£  -  30 


(4.1) 


which  idealizes  the  result  of  embedding  a  trace  of  molecular  oxygen  in  a 
bath  of  hot  atomic  oxygen  (the  dominant  reaction  is  dissociation). 

The  rate  equations  are 


dt^  =  'Kf  "o  V 


(4.2) 


+2Kf  no  V 


(4.3) 


which  have  the  conservation  law 


n„  +  2n  ~  =  c  =  const.  >  0 
o  oZ 


(4.4) 


corresponding  to  the  conservation  of  the  total  number  of  oxygen  atoms,  free 
and  bound  in  pairs.  Equation  (4.3)  can  be  rewritten  as 


ar '  2Kf\V  =  Kfno(c '  "o’ 


(4.5) 


Introduce  the  number  fractions 


Then,  (4.4)  implies 


B -f  •  ’'‘Hr 


A  +  B  =  1 


and,  since  A  and  B  are  positive. 


(4.6) 


(4.7) 


0<A<1  ,  0  <  B  <  1 


(4.8) 
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We  rescale  for  simplicity  the  time  by  K^c  and  obtain  from  (4.5) 


^  =  -AB  =  -A  +  A2  (4.9) 


The  last  equality  follows  from  (4.7).  From  (4.9)  we  find  the  equivalent 
pa i rs  : 


and 


j  =  -I  +  7 


7  =  -2A2  +  2A3 


-  -  T 

B  =  +B  -  B 


“T  ~7  “T 
B  =  2B*  -  2BJ 


(4.10a) 

(4.10b) 

(4.11a) 

(4.11b) 


First-order  closure  consists  in  using  (4.10a)  together  with  a  model 

A2  =  F(A) 

"T 

We  apply  the  bounds  on  A  derived  in  the  text  (Section  II) 


(4.12) 


-2  77- 

A  <  A  <  A 


(4.13) 


and  prove  the  following. 

(I)  Given  A(0)  ,  the  solution  of  any  model 

t  =  -A  +  F(A)  (4.14) 

where  F  satisfies  the  bounds  of  Eq.  (4.13),  is  bounded  above  by  the 
solution  of 

t  =  -A  +  A  =  0  (4.15) 

and  bounded  below  by  the  solution  of 

j  =  -A  +  A2  (4.16) 

(II)  Equation  (4.15)  gives  the  exact  solution  of  (4.9)  with  the  initial 
probability  density 
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1 


P(Aq)  =  (1  -  Ao)6(Aq)  +  Aq6(A0  -  1)  (4.17) 

and  Eq.  (4.16)  gives  the  exact  solution  of  (4.9)  with  initial  density 


P(A  )  =  6(A  -  A  ) 
'o'  v  o  o 


(4.18) 

Subscript  is  used  to  denote  the  time  at  which  the  variable  is  considered. 

(Ill)  There  exists  a  nonextremal  choice  of  F  which  does  not  yield  the 
exact  solution  of  (4.9)  for  the  corresponding  density.  The  exact  solution 
of  (4.9)  is 

.-t 

(4.19) 


At  = 


A  e 
o 


1  "  +  Art  e 

o  o 


-t  -  *'1<Ao> 


The  initial  value  of  A,  Aq  is  conveniently  expressed  in  terms  of  the 
value  of  A  at  t  by  inverting  (4.19) 

+t 


A0  =  - — - Tf  -  *(At) 

1  "  At  +  At  e 


(4.20) 


which  can  also  be  obtained  by  time  reversal  of  (4.9). 

We  now  consider  the  above  points  successively.  We  write  the  general 
solution  of  (4.14)  as 

e-(t-X)  dx 

0 


(4.21) 


where  A  is  understood  to  be  a  function  of  A  .  Since 

e-(t-A)  > q  a-|i  A 


we  have  from  (4.13U) 


and  from  (4.13L) 


/: 

r 


,-(t-A) 


a-(t-x) 


7(X)  -  A(X)  dX  <  0 


?(X)  -  A2(X)  dX  > 


(4.22) 


(4.23) 


(4.24) 
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Hence, 


* 

f  e'^t"X^A2(X)  dA  <  f  e'^t'X^7(\)  d A  <  f  e‘^'A^A(A) 

Jn  Jn  Jn 


dA  (4.25) 


-t; 


If  we  add  e  Aq  to  each  term  of  the  double  inequality,  we  find 

^  t 

~  (  t— A )  «  2  y 


e-tA0  +  f  e'^'^A1 

•'o 

ie'‘So  *f  1 

Jo 

*J  < 

JQ 


(A)  dA  <_ 
-^-X)7(A)  dA  < 
■(t~X)A{A)  dA 


(4.26) 


which  proves  our  first  proposition;  i.e.,  the  solutions  of  the  extreme 
models  are  the  bounds  of  the  solution  of  any  (statistically  allowed)  model. 
Furthermore,  the  exact  solution  (4.19)  when  averaged  with  (4.17)  gives 


At  =  Ao 


(4.27) 


which  is  indeed  the  solution  of  (4.15).  Also,  (4.19)  when  averaged  with 
(4.18)  gives 


A.  e 


-t 


At  = 


1  -  A„  +  A„  e 
o  o 


(4.28) 


which  can  be  readily  verified  to  be  the  solution  of  (4.16).  We  have  thus 
proven  our  proposition  II  and  thus  demonstrated  that  the  extreme  models 
represent  exact  solutions.  Therefore,  combining  I  and  II,  we  conclude  that 
any  statistically  correct  choice  for  the  model  F(A)  will  remain  statis¬ 
tically  correct. 

A  simple  example  of  III  is  the  following.  Model  F(A)  by  the  mean 


7 


Mod  1-?  i  - 
2  A  +  2 


(4.29) 


A  corresponding  density  is 
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P<V  •  ?  «(Ao  •  V  *  i 


(1-A0)5(A0)  +  A0«(A0-  1) 


(4.30) 


A  short  calculation  shows  that  the  solution  of  the  model  equation 


-a  + 


~  A^  +  A 


is 


Ao  e 


-t/2 


i  -  a  +TT177 

O  0 

while  the  average  of  (4.19)  with  P(A  )  given  by  (4.30)  is 


f  A  e-t 

At  =  I  - 2 - It 

J0 


Ao  r 

P(Ao)  dAo  - 


-t 


1  + 


1  -  A  +  A  e 
o  o 


-t 


(4.31) 


(4.32) 


(4.33) 


which  does  not  coincide  with  (4.32)  although  (4.32)  is  properly  bounded  as 
discussed  above. 

The  bounds  obtained  on  At  from  first-order  closure  (i.e.,  given  AQ 
only)  are  loose  bounds.  We  can,  however,  proceed  to  second-order  closure 

3 

and  tighten  them  by  giving  AQ  and  AQ  .  The  statistical  bounds  on  A 
(from  Section  II)  are 


_ 

n 


2, 2 


AL.  <  a3  <  A2  -  (A  -  A_) 
A  1  -  A 

1? 


(4.34) 


The  equations  for  A  and  tc  are  given  in  (4.10a)_and  (4.10b).  Second- 
order  closure  corresponds  to  choosing  a  model  for  A 


A3  =  G(A  ,7) 


(4.35) 


which  we  assume  subject  to  the  bounds  (4.34).  We  choose  G  to  correspond 
to  the  lower  and  upper  bounds  in  (4.34).  For  the  lower  bound,  we  find  that 
the  quantity 

72 

(4.36) 


A2 

KL  =  = 

•  2 
PC 


is  a  constant  of  the  motion.  Integration  then  yields  (AQ  and  Aq  given) 
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(4.37) 


A2  e"* 


At  = 


~~7  ~7  -t 

a  -  A  +  A  e  x 
0  0  0 


The  density  corresponding  to  (4.34L)  is  (from  Section  II) 


L  •  2 


P(A0)  -rf  *(V*4<Ao-f 

Ao  A  '  Ho 

0 


(4.38) 


A  simple  calculation  shows  that  the  exact  (4.19)  averaged  with  (4.38)  gives 
(4.37). 

The  upper  bound  is  treated  conveniently  in  terms  of  4.11a)  and  (4.11b). 
In  this  case,  the  constant  of  the  motion  is 


„  _  B2  _  1  -  2A  +  A2 

u  T  ~1 

a  1  -  2A  +  K 


(4.39) 


The  corresponding  density  is 


P(A0>  -  ^  «(Aq  -  1)  A  A  -  ^ 

K  K,  x  ° 1 


(4.40) 


and  the  model  solution  coincides  with  the  average  of  (4.19)  over  (4.40)  as 

(4.41) 


B0  e+t 


Bo  -  %  *  Bi  e+t 


which  can  be  written  as 


At  = 


A2  -  A2  +  (A  -  7)  e_t 
0  0  0  0 


1  -  +  ^  +  (A„  -  7)  e_t 

0  0  0  0 


(4.42) 


The  bounds  given  by  (4.37)  and  (4.42)  are  substantially  more  stringent  than 
those  given  by  first-order  closure. 

As  an  example  of  nonextremal  choice  of  second-order  closure,  we  consider 
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(4.43) 


A3  =  (1  +  A)(/?  -  A2)  +  A3 


corresponding  to  the  density  (frozen  A  model  of  Section  II) 


A1  ^  a“ 

P<A  )  ■  ~~  6{A  -  1)  t  . 

Bo  \ 


A' 


AoBo 


6(Ao  -  Aq) 


(4.44) 


The  solution  of  the  model  equation  is  obtained  by  exploiting  the  integral 
resulting  from 

t -  - 25  <4-45> 

and  is 


t  "  2 


1  T  -2 
i.  .  A  +  A  x 
4oo 


An  -  i  - 
o  2  * 


I  _  a2  +  A2  +  e‘2t 
4  o  o  e 


A0  -  I  V  T  “  ^  +  Ao 


An  -  A-  - 
o  2  » 


1  «2  .  7  2  -2t  I  -n 

4  "  o  Ao  '  e 


while  the  exact  solution  (4.19)  averaged  over  (4.44)  yields 


A’. 


At  = 


Ao  e 


-t 


AoBo/ 


1  -  A  +  A  e 
o  o 


-t 


(4.46) 


(4.47) 


< 
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Figure  5.  Interpretation  of  §  axis  on  a  stretched  y/(x-x  )  axis. 
This  interpretation  eliminates  y(x  -  x  ).  0 
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Figure  3.  Three  important  special  cases  of  Figure  7 
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Figure  24.  Given  mole  fraction  moments,  obtain  bounds  on  mass  fractions 
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Figure  26.  Probability  density  function;  r  =  0.38 


Figure  27.  Failure  of  a  square  hat  to  fit  a  triangular  distribution 
(finite  probability  of  negative  molar  fraction) 


Figure  28.  Allowed  domain  for  symmetric  models 
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Figure  30.  Probability  density  function;  r  =  0.38, 
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Figure  33.  Probability  densi 
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Figure  37.  Normalized  third  moment  compared  to  experiment  for  A 
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Figure  43.  Normalized  third 
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Figure  44.  Predictions  for  pure  Ng  compared  to  experiment 
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Figure  45.  Predictions  for  pure  He  compared  to  experiment 


iction;  r  =  0.38;  s  =  7.0.  Model  with  four  cells  and  largest  entropy 
poorness  in  disagrement  with  experiment. 
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Figure  49. 
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Third  moment  compared  to  experiment  for  pg/pj  =  7.0,  u2/u^  =  0.38 
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